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^ ABSTRACT 
O 

The clustering in redshift space is studied here to first order within the framework of 
i-H gravitational instability. The distortion introduced by the peculiar velocities of galaxies 

C*~) results in anisotropy in the galaxy distribution and mode-mode coupling when analyzed in 

Fourier space. An exact linear calculation of the full covariance matrix in both the real and 
^ Fourier space is presented here. The explicit dependence on fio and the biasing parameter 

On 

is calculated and its potential use as a probe of these parameters is analyzed. It is shown 

Oh 

^ that Kaiser's formalism can be applied only to a data set that subtends a small solid angle 

S_i 

^ on the sky, and therefore cannot be used in the case of all sky surveys. The covariance 

a 

matrix in the real space is calculated explicitly for CDM model, where the behavior along 
and perpendicular to the line of sight is shown. 
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I. INTRODUCTION 

One of the basic tools of observational cosmology is the mapping of the galaxy distribu- 
tion in redshift space, where the observed redshift is used rather than physical distances. 
Indeed, redshift surveys have played a major role in shaping our understanding of the 
large scale structure of the universe. Given the tremendous difficulties in obtaining true 
distances of galaxies, redshift surveys will remain the major method of studying the large 
scale distribution of galaxies. Various groups have determined the statistical properties 
of the galaxy distribution directly from galaxy redshift catalogs. In particular, the power 
spectrum has been evaluated from the galaxy redshift distribution (Boumgart and Fry 
1991, Park, Gott and da Costa 1992, Fisher et al. 1993, Feldman, Kaiser and Peacock, 
1993). 

The most basic analysis of the large scale clustering of galaxies is that of the power 
spectrum estimation, or equivalently its Fourier conjugate the 2-point correlation function 
(cf. Peebles 1980). From the theoretical point of view the power spectrum can be easily 
calculated in a variety of models and thus the observationally estimated power spectrum 
can be used to yield significant information on the nature of the early universe and the 
process of structure formation. However, theoretical calculations usually give the real 
space power spectrum while observationally it is determined in redshift space. Now, the 
formation of structure out of an expanding Friedmann universe involves departure from 
homogeneity, and therefore by the continuity equation peculiar velocities arise. This leads 
to the displacement of galaxies along the line of sight and introduces anisotropy to the 
apparent clustering of galaxies. On small scales, corresponding to nonlinear perturbations, 
virial velocities in collapsed clusters leads to the so-called 'finger of god' effect (Peebles 
1980). This results in the elongation of the structure along the line of sight in redshift 
space. On larger scales where deviations from the Friedmann universe are in the linear 
regime, the infall to overdense regions compresses the structures along the line of sight. 
The relation between the power spectrum in real and redshift space is model dependent 
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and can be solved only within the framework of an assumed model. 

At the present paper we shall focus on the theoretical calculation of the redshift space 
power spectrum in the linear regime. This was first calculated by Kaiser (1987) who 
used the single mode approximation to find a simple relation between the redshift and 
real space power spectra. This relation depends on the density parameter fio and the 
biasing parameter b (see discussion below), and thus can be used to find these important 
parameters. Kaiser's simple relation is rigorously correct in the case of a single Fourier 
mode, but as is shown in the present paper it is not valid in the general case of a multi- 
mode field which subtends a wide solid angle over the observer's sky. Nevertheless, Kaiser's 
formalism has been widely used by many workers in the field to reconstruct the real space 
power spectrum from the observed one and to set constraints on fio and b. Therefore it 
cannot be applied to analyze the (almost) full sky redshift surveys such as the IRAS 1.2 Jy 
redshift survey (Fisher et al. 1993) and the IRAS QDOT redshift survey (Feldman, Kaiser 
and Peacock 1993). This limitation has escaped the attention of many workers in the field 
and has been used in analyzing full sky surveys. The major aim of the present paper is to 
clarify this important point and to present a rigorous calculation of the correlation function 
and the power spectrum in redshift space. Here we shall focus only on the theoretical 
aspects of this calculation and in a forthcoming paper this will be applied to the IRAS 
data base (Bistolas, Zaroubi and Hoffman, 1994). The exact treatment of selection effects 
and shot noise will be given at that paper. The main results reported here are, first, that 
the anisotropy introduced by the redshift manifests itself by an off-diagonal coupling of the 
Fourier modes, where the Fourier space is the conjugate to the redshift space and the ratio 
of the power spectra of the two spaces is ^-dependent. Thus, the power spectrum is not the 
Fourier transform of the two-point correlation function. Here, the term power spectrum is 
used in the more general sense of the covariance matrix in the Fourier space. The second 
result is that the correlation between any two points in the redshift configuration space 
depends not only on the magnitude of the displacement vector which connects them but 
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also on its orientation and distance from the observer. 

The structure of the paper is as follows. An exact (to first order) calculation of 
the Fourier modes covariance matrix in redshift space is presented in §11. The Kaiser's 
formalism is recovered here in the special case where the data subtends a small solid angle 
on the sky (§111). The covariance matrix in the redshift configuration space is presented 
and calculated explicitly for special cases in §IV. A general discussion is given in §V. 

II. LINEAR THEORY IN REDSHIFT SPACE 

The coordinate transformation from r-space (real space) to s-space (redshift space) is 
given by: 

U(r) - U(0), 

s r = r 1 + ^ 1 

r 

Here U(r) and U{0) are the projection of the peculiar velocity at the point r and at the 

position of the observer, respectively, on the radius vector connecting these two points. 

Distances are measured here in velocity units. The relation between the densities in the 

two spaces are given through the Jacobian, 

d 3 s U(r)-U(0) dU(r) 
Ps = Pr| f = ^ 1 + ~ ] + • (2) 

The derivation here assumes that the selection effects are given by the selection function 
(f)(r) and the measured density in redshift space is weighted with (/> _1 (,s). To first order (in 
perturbation theory) the fractional overdensities in s- and r-space are related by: 

Wr) = 4r(r) _( 2 + ^)EW^_^ (3) 

\ a(lnr) J r or 

Note that in the second term of Eq. 3 (U(r) — U(0)^j is of the order of the rms peculiar 
velocities and it is distance independent. This second term scales as the inverse of the 
depth of the survey and is usually neglected as compared to the other terms (Kaiser 1987). 
However, in our patch of the universe velocities are of the order of several 10 3 km/ sec which 
is not negligible on scales of ~ 10 4 km/ sec. Therefore, our treatment will in general include 
this term unless it is stated otherwise. 



The relevant Fourier transformed fields are: 



8 e (r) = J exp(ir -k)^(k)d d A; (4a) 

8 r (r) = y exp(ir • k)5 r (k)d 3 £; (46) 
From linear theory (Peebles, 1980) one finds 

U = -if (tl ) [ exp(ir • k) ^ S r (k)d 3 k 
J kr 

BTT f k • r 

— = -f(n ) J exp(ir • k )( — f Sr (k)d 3 k, (5) 

where /(f^o) ~ ^o' 6 - I n the case of linear biasing /(f^o) is to be replaced by (3 = f(Q,o)/b. 
In the inverse &-space Eq. 3 is written as: 



6 8 (k) = s r (k) + vi (k) + f(n ) J d s k's r (k')(i(k',k)- tl2i ^ k,) ) 



(6) 



(7c) 



where io(k) and the coupling kernels J(k', k), /i(k', k) and l2(k',k) are defined as follows: 

J (k) = /exp(-tk -r)(2 + (7a) 

K ' (2tt) 3 J y V d{lnr) J kr y ' 

/(k '' k)= (^ / d3 ^Mir-(k'-k))(^ f (76) 
J 1 (k',k) = ^(k-k') = ^ ? | exp^k'-k)-!-)^ 

^ k) = / exp(i(k ' - k) • r) ( 2 + W) ) ^ (w) 

The crucial new result obtained here is that the transformation to redshift space 
variable, 6 S , introduces mode-mode coupling . Now, in the standard model of cosmology 
the perturbation field is assumed to be Gaussian, and for a homogeneous and isotropic 
field, this implies that the covariance matrix of the real space density field is diagonal, 

<8 r (k)8* r (k')>= P(k)8 D (k-k') (8) 



and the covariance matrix is given by: 



<^(k 1 )^:(k 2 )> 



= J <z 3 jfe;p(k'i) 




/ 2 (k' 1 ,k 1 ) + /(Oo)/(k' 1 ,k 1 )] 



(9) 



x [J 1 *(k' 1 ,k 2 ) + 



/ 2 *(k' 1 ,k 2 ) + /(n )r(k' 1 ,k 2 )] 



+ F 2 / (k 1 )/*(k 2 ) 



Neglecting the second term in equation (3) the covariance matrix thus obtained is, 



<S e (k)S* e (k')> = P(k)S D (k-k') 



+ f(n )p(k')i(k',k) + f(n )p(k)r(k,k') 

+ f(n ) [ d 3 A;"P(A ; ")/(k",k)r(k",k') 



(10) 



which still includes a mode-mode coupling and it is much more complicated than the simple 
formula commonly used. 

The transformation to redshift space changes the real space perturbation field in two 
fundamental ways. The first and obvious one is the anisotropy introduced by the line- 
of-sight peculiar velocity. This was recognized by Kaiser (1987) and in all subsequent 
works. The second and more subtle effect is the mode-mode coupling in the Fourier space 
representation. The source of this coupling is the geometric nature of the real - redshift 
space transformation which depends on the term. Many of the common properties of 
the real space perturbation field do not hold in redshift space. One such basic relation 
often used is that the real space 2-point correlation function and the power spectrum form 
a Fourier conjugate pair. This relation breaks down in redshift space as is manifested by 
Eq. 9 which shows that the statistical properties of the covariance matrix (i.e. the power 
spectrum) only. In any case, there is no simple ^-independent relation between the real and 
redshift space power spectra. This is clearly manifested by the work of Fisher et al. (1993), 
which generated mock IRAS catalogs from cold dark matter model (CDM; Blumenthal 
et al. 1984) N-body simulations and compared the power spectra evaluated in the two 
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spaces. Figs. 3b and 3c of Fisher et al. clearly shows a k dependence of this relation, and 
in particular that this is inconsistent with Kaiser's relation. This result was also confirmed 
by Gramann, Cen and Bahcall (1993), who have analyzed N-body simulations in redshift 
and real space. They calculated the power spectrum in both spaces and found that the 
ratio of these is k dependent. This is manifested by Fig. 2a of Gramann, Cen and Bahcall 
(1993). Last point to be made here concerns a property which holds also in redshift space. 
The redshift transformation is linear and therefore it preserves the Gaussian nature of the 
random field. Now, the 6 s (s) field is not isotropic and its (&-space) covariance matrix is 
not diagonal but its statistical properties can be analyzed within the standard framework 
of Gaussian random fields. 



The single mode approach of Kaiser (1987) is easily recovered from our more general 
formalism, we apply his assumptions to equation (6), namely by neglecting lo and I2, and 
substituting 



III. COMPARISON WITH KAISER'S FORMALISM 



S r (k) = S S D (k - k ) 



(11) 



we obtain Kaiser's result, 




(12) 



Expressing Eq. 12 in Kaiser's notation we define 



k • r 



(13) 



k r 



to obtain: 



s 8 (k ) = s (ko)(i + f(n )fi 2 ) 



(14) 



Here the angle- averaged of ji 2 



is ji 2 = 1/3. Kaiser's basic relation is 



tf a (k ) = tfo(ko)(l + /(fto)/i 2 ), 



(15) 
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without taking the angle- average. However, this expression is mathematically ill defined, 
as any r dependent quantity does not exist in k space. Only after the angle- average is 
taken, namely integrating over real space, Kaiser's results is properly defined and is equal 
to the expression obtained here. In this sense Kaiser's result has been used by workers in 
the field (eg. Fisher et al. 1993, Feldman, Kaiser and Peacock, 1993), without realizing 
that it is rigorously correct only in the case of a single mode. 

Kaiser's formalism is useful and accurate in the case where the data is confined within 
a volume which extends a small solid angle relative to the observer. In this case the //. is 
redefined to be: 

ko • (R + r) 
k Vr 2 + R 2 

Here R defines the center of that volume, in some general sense. Now, in the above limit 
of small solid angle, R ^> r, and 

"~iwr (17) 

Substituting this approximation of //. in Eqs. 5-10, the mode-mode coupling that is intro- 
duced by //. breaks down and Kaiser's results are fully recovered, namely Eq. 15 and 

P s (k) = P(k)(l + f(Sl )v 2 ) 2 . (18) 

Indeed, this is the case studied by Cole, Fisher and Weinberg (1993), who analyzed redshift 
distortions in N-body simulations and confirmed Kaiser's formalism in the limit of samples 
that subtends a small solid angle. 

The exact correlation matrix (Eq. 10) has a complicated mode-mode coupling, which 
is introduced by the J(k,k') coupling term. It can be shown that for the diagonal terms 
in the limit of small wave number (k = k' — > 0) one obtains: 

<|^(k)| 2 >~P(k)^l + ^/(O )j (19) 

Note that this differs from Kaiser's formula in not having the quadratic /(fio) 2 term, which 
for /(f^o) = 1-0 amounts for a ~ 12% difference. It is interesting that the asymptotic 



relation found here is the same result quoted by Kaiser as the one derived from Peebles' 
(1980) formalism. 

IV. CORRELATION FUNCTION IN REDSHIFT SPACE 

As discussed before the transformation to redshift space introduces anisotropy and 
mode-mode coupling in Fourier space representation. This complicates the picture in 
Fourier space and motivates us to calculate the correlation matrix in the redshift configu- 
ration space. The correlation matrix is expected to have three degrees of freedom (d.o.f), 
the first is the distance of one of the points from the observer and the other two are, the 
line of sight (l.o.s.) component of the displacement vector which connect the two points 
and the component perpendicular to it. The correlation function in redshift space for the 
points i"i and r 2 can be calculated directly from Eq. 3. Let r = r2 — i"i, r = |r| and r = ^, 
the correlation matrix is given by, 



T" ( v 1 T" ( v 1 

< £ s ( s i) S s (s 2 ) >= £(r) — r 2 • »?(r) H — ri • f?(r) - r 2 • dri(r)/dr 2 

r 2 n 

+ ri • dr, r ) d n + V 1J V 2} r}tf(r)r 2 + r\— ^r 2 

ri r 2 ri dr 2 

T(r 2 ) 4 M(r)„ 4 9 2 4(r)„ T(n) T(r 2 ) , „ 
+ ti^T 2 ^ + r + V U V 2 > V 2 co S 8 1 co S 8 2 20 

r 2 or i or\Or 2 r\ r 2 

where T(r) = 2 + dln((f)(r)) / dln(r) . 

The tensor $ = {^ij} is the two-point correlation tensor of the velocity field (Gorski 
1988), 

< vi(n) vj(r 2 ) >= tf^r) = *±(r)8 i>j + [tf „(r) - *±(r)]r i r j (21) 

where ^(t") and \t ±(r) are the radial and transverse velocity correlation functions respec- 
tively. The spectral representation of these correlation is given by 

*-L.ll( r ) = H2 ° f 2 f 2 ° )2 jT P(k)K ul (kr)dk, (22) 
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where K±(x) = ji{x)/x and K\\(x) = jo — 2j\(x)/x and ji(x) is the spherical Bessel 
function of order I. 

The vector rj = {r/i} is the two-point cross-correlation vector between the density and 
velocity fields, and it is given by 

< Yi(rO *(r 2 ) >^ ifc(r) = g °^" o) B(r)r,- (23) 

where H(r) = J °° kP(k)j\(kr)dk. Note that the vector 77 is parallel to r which is the 
natural direction in this case. 

Numerical calculation of Eq. 20 is rather straightforward. Here we present calculations 
of this equation for an 8Mpc/h Gaussian smoothed CDM model with no selection effect 
and with Vo = 0. Figs, l(a-c) show plots of the correlation function in redshift space 
around a point at a distance of 55, 100 and 500 Mpc/h from a given observer, respectively. 
Fig. Id is the correlation function calculated by Kaiser's (1987) formula, (1 + 2/3/(fio) + 
l/5/ 2 (fio)) C( r )- The Z axis is taken to be the l.o.s. direction. Figs, l(a-c) demonstrate 
the different behavior of the correlation function, parallel and perpendicular to l.o.s.. They 
also show the effect of the distance from the observer on the correlation function, which 
is caused by the U/r term which indicates that it can't be neglected in the analysis of the 
nearby Universe 100 Mpc/h). 

Fig. 2 shows a linear-log plot of the correlation function in the direction of the 
l.o.s. and perpendicular to it, calculated from Eq. 20 and compared with the real-space 
correlation function and with Kaiser's (1987) approximation (angle averaged). The two 
figures show that the zero crossing of the correlation function in redshift space occurs at a 
distances shorter than the ones expected from the real space correlation function or from 
Kaiser's approximation. This can be simply explained by the mass conservation, admitted 
under the redshift transformation, which amplifies the clustering as r — > and therefore 
attenuates it as r 2 recedes from ri. 

This result is confirmed by the the work of Lilje and Efstathiou (1989) where they used 
standard CDM simulations to calculate the correlation function along the l.o.s. for galaxies 
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within given transverse separation range (in particular their Figs. 4-7). Fig. 4 of Lilje and 
Efstathiou shows the N-body calculated correlation function in the real and redshift space 
along the l.o.s. for galaxies with transverse separation in the range of 0.5-1.3 Mpc/h. Note 
that they ignored the variation of the redshift space correlation function with the distance 
from the observer. Their Fig. 4 confirms our findings with regard to the l.o.s. correlations 
shown in here in Figs. 1 & 2. Namely the l.o.s. redshift-space correlation function is larger 
than the l.o.s. real-space correlation function on small scales and is smaller than it on large 
scales, where the transition occurs on the scale of 10-20 Mpc/h (for standard CDM). This 
behavior repeats itself in their Figs. 5-7 where the l.o.s. redshift-space correlation function 
is calculated for a larger transverse separation. Lilje and Efstathiou results also agree with 
the general behavior found here that the zero crossing scale of the redshift-space function 
is smaller than that of the real space covariance function. It also disagree with Kaiser's 
approximation where the ratio of the ratio of the two functions is constant (> 1). 



V. DISCUSSION 

A proper, first order, calculation of the relation between the actual real space density 
field correlation function and the one observed and measured in redshift space has long 
be overdue. The formal tools of the the linear theory of gravitational instability could 
have been easily applied to improve on Kaiser's (1987) approximate results. Given the 
importance of the redshift-to-real space transformation it is rather surprising that the 
problem has not been yet properly addressed. The present paper fills up this gap and 
introduces a straight forward linear calculation of the redshift space covariance functions 
in the configuration and Fourier spaces. 

The present calculation serves a dual purpose. On the one hand this can be used to 
develop an efficient algorithm for the reconstruction of the real space density field from 
the observed redshift galaxy distribution. This can be done in the framework of Wiener 
filtering (Zaroubi et al. 1994, Fisher et al. 1994 ). On the other hand the redshift space 
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covariance matrix can be used to estimate relevant cosmological parameters by maximum 
likelihood analysis (eg. the power spectrum, fio and the biasing parameter). Now, the 
statistical properties of the redshift space perturbation field differ from real space field in 
two major ways, namely the transformation induces mode-mode coupling and anisotropy 
of the covariance matrix. The isotropy and orthogonality of the linear regime of the k- 
space covariance matrix is expected in a wide range of cosmological models, and therefore 
redshift space deviations from these can be used as a probe on fio or (3, in the case of 
linear bias. However, in order to make the present calculations into a useful cosmological 
probe one should complement these by numerical N-body simulations to study the role of 
non-linear effects and to establish the dynamical range over which the linear theory that 
has been developed here applies ( Bistolas, Zaroubi and Hoffman, work in progress). 
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FIGURE CAPTIONS 

Fig. 1: A) The correlation function for the a standard CDM model with 8Mpc/h smoothing, 
as measured in the redshift space according to the exact linear calculation around point 
at 55 Mpc/h distance. Here we take <f)(r) = 1 and Vo = 0. All contours are spaced at 
0.05 with solid (dashed) lines denoting positive (negative) contours. The heavy solid 
line denotes the zero line. B) The same as A for point at 100 Mpc/h distance. C) 
The same as A for point at 500 Mpc/h distance. D) The same correlation function as 
calculated by Kaiser's (1987) approximation; £ s = (1 + 2/3/ + l/5/ 2 )^ r . 

Fig. 2: The correlation function of the same model described in Fig. 1, parallel to the l.o.s. 

(dotted-dashed line) and perpendicular to it (solid line). The correlation function 
according to Kaiser's (1987) approximation is also shown (dotted line) and the real 
space function is given by the dashed line. Make sure the lines are marked 
correctly, after my correction! 
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Figure 1 
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